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Abstract. Increasing evidences provided primarily by 
the cluster lensing and numerical simulations of cluster 
formation indicate that galaxy clusters may contain com- 
pact cores that are substantially smaller than the cores 
revealed by the X-ray observations of hot intracluster gas. 
In this paper we present a model that describes how two 
distinct cores can grow simultaneously as a result of in- 
fall from the background dark matter. This model needs a 
pre-existing localized large fluctuation, which can be the 
non-Gaussian density peak of the primordial fluctuations, 
to seed the individual growing cluster. On the other hand, 
numerous recent observations also show that no strong 
evolution is detected for galaxy clusters within the red- 
shift up to 2 ~ 0.8. We therefore present a comparison of 
observations with the saturated cluster configuration re- 
sulting from our infall model. In the saturated state, the 
predicted compact core mass is about fewxlO^^A/0 and 
the core size about 1 kpc for a cluster mass ~ several 
X 10^^ Mq within a radius of 3 Mpc. We have successfully 
reproduced the dark matter distribution revealed by the 
gravitational lensing, and the observed radial distributions 
of cluster galaxies, intracluster gas (i.e., the conventional 
P model) and baryon fractions in a consistent way. This 
model, when combined with the observed X-ray surface 
brightness profiles of clusters, predicts that the overall 
temperature of intracluster gas has a tendency of radial 
decline with a mean polytropic index 7 « 1.2. Finally, 
the so-called f3 discrepancy finds a natural explanation 
for such a polytropic gas in our model. 

Key words: cosmology: theory - dark matter - galaxies: 
clustering — hydrodynamics 



1. Introduction 

The existence of a compact and massive dark-matter core 
in the galaxy cluster, which is much smaller than the X-ray 
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gas core of the empirical /3 model, has recently acquired 
substantial supports from various evidences. It has been 
recognized since the early study of the very thin arclike 
images of distant galaxies created by the potentials of the 
foreground lensing clusters that the dark matter profiles 
of the lensing clusters should be sharply peaked towards 
cluster centers (e.g. Hammer 1991). In order to account 
for the total number of observed giant luminous arcs, the 
core size of dark matter distribution modeled by a softened 
isothermal sphere should be at least ~ 10 times smaller 
than the one for luminous matter profile (Wu & Hammer 
1993). Grossman & Saha (1994) strengthened this argu- 
ment by explicitly setting an unprecedentedly small upper 
limit on cluster core radius, rdark,c < 11 kpc, for a soft- 
ened isothermal sphere. Many recent work on modeling of 
arcs/arclets have essentially confirmed these findings (e.g. 
Giraud 1999; Wilham, Navarro & Bartelmann 1999). Even 
for the X-ray measurements of hot intracluster gas, the dy- 
namical analysis also arrives at a similar conclusion that 
the binding matter distribution of clusters should have a 
core radius of rdark,c < 100 kpc (Durret et al. 1994). Op- 
tically, the omni-presence of cD galaxies in the cores of 
many rich clusters may indicate unusual accumulation of 
central masses. Finally, a singular mass profile of p cx r^^ 
at cluster center is predicted by the numerical studies of 
cluster formation in the framework of the standard CDM 
model as well as other cosmological models (e.g. Navarro 
et al. 1995; Cole & Lacey 1996; Navarro, Frenk & White 
1997; Fukushige & Makino 1997; Eke, Navarro & Frenk 
1998). Such a singular mass profile also houses a weak 
compact core. 

That is, while the cluster X-ray gas often shows a core 
of typical size about 100 — 250 kpc, the dark matter core 
can instead be substantially smaller than 100 kpc. The 
co-existence of two distinct mass cores in a single clus- 
ter appears to have challenged our current understanding 
of cluster formation. Although a great deal of theoreti- 
cal efforts have in the past been put forth attempting to 
grasp an understanding of the cluster physics, few have 
considered the possibility of a massive compact dark mat- 
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ter core in the galaxy cluster. Past analytic investigations 
for the gravitational infall of matter into galaxy clusters 
have been pursued mainly toward the following two is- 
sues: the cosmic evolution of cluster properties such as 
the total mass, number counts, gas luminosity and tem- 
perature (e.g. Gunn & Gott 1972; Press & Schechter 1974; 
Kaiser 1986; 1991), as well as the gas and dark matter 
density profiles inside clusters (e.g. Fillmore & Goldre- 
ich 1984; Bertschinger 1985; Teyssier, Chieze, & Alimi 
1997). Regardless of their oversimplification, these ana- 
lytic treatments based on self-similar models indeed pro- 
vide the gross properties of clusters that were shown to 
be in gross consistency with optical/X-ray observations 
(Kaiser 1991). Moreover, these analytical studies can also 
reveal the underlying cluster evolution mechanisms, in a 
complementary manner to those shown by the particle 
and hydrodynamic simulations. In spite of the past suc- 
cess, previous studies commonly obtained (or assumed) a 
single core within the cluster. 

In the self-similar regime, clusters are, however, shown 
to undergo very rapid evolution (Gunn & Gott 1972; 
Kaiser 1986). When these results arc combined with the 
recent observations made primarily at X-ray wavebands, 
which demonstrate that no evolution has been detected 
within the redshift at least out to z « 0.8 among the clus- 
ter number counts. X-ray luminosity, temperature func- 
tion as well as typical scale lengths (e.g. Fan, Bahcall & 
Cen 1997; Rosati et al. 1998; Vikhhnin et al. 1998; Wu, 
Xue & Fang 1999 and references therein), one is led to a 
conclusion that the evolution of these clusters must have 
ended at an early epoch. 

In this paper, we will concentrate on the dark mat- 
ter and the galaxy/gas distributions of clusters revealed 
mainly by the saturated state of a new class of self-similar 
infall model, which describes the simultaneous build-up of 
two distinct cores in a single cluster. Our starting point 
differs from the past theoretical models in the following 
way: (a) Unlike the previous work which chose a power- 
law initial perturbation in the mass (e.g. Fillmore & Gol- 
dreich 1984; Teyssier et al. 1997), we seek a solution to the 
infall problem of dark matter particles by assuming an ini- 
tial compact seed density fluctuation at the site of grow- 
ing cluster. This aspect is in line with the works of Gott 
(1975) and Bertschinger (1985), and the large-amplitude 
seed fluctuation may arise from a non-Gaussian peak in 
the primordial density fluctuations, (b) Instead of assum- 
ing the collisionless particles to move with radial orbits 
(Merritt, Tremaine & Johnstone 1989), which are known 
to be unstable (Binncy & Tremaine 1987), we consider 
the other extreme situation where the particle orbits in- 
side the cluster are so chaotic that the velocity dispersion 
is essentially isotropic. 

After a substantial growth of both cores, they can 
rapidly reach a saturated state, in which the mass ratio of 
the two cores remains approximately the same as that dur- 
ing the infall phase. (This is because both cores are gravi- 



tationally bound, and they have been grossly virializcd re- 
gardless of whether there exists an infall.) Our attempt at 
comparing the theoretical predictions with observational 
results will therefore be made at this late stage when the 
matter infall has ended, and a hydrostatic equilibrium has 
already been established between the dark-matter gravi- 
tational potential and the cluster galaxies and intracluster 
gas. 

We begin in section 2 by presenting a specific model 
describing the self-similar evolution of spherical collapse of 
collisionless cold dark matter in an open and a closed uni- 
verse. (The former also represents the situation where the 
infall began only recently and is still ongoing in a flat uni- 
verse.) We then derive the density profiles of dark matter 
in the infall phase and in the saturated phase, respectively. 
In section 3 the predicted dark matter distribution and the 
galaxy /gas density profiles are compared with those re- 
vealed by the gravitational lensing and optical/X-ray ob- 
servations. Finally, the discussion and summary arc given 
in section 4. Beginning from section 3, where comparisons 
of the theoretical predictions with the observations are 
made, we work mainly with the saturated state (except 
in section 3.1), and the Hubble constant is chosen to be 

= 50 km s~^ Mpc~^ for convenience. Throughout this 
work, we set the cosmological constant equal to zero. 

2. Self-similar solutions for the collapse of cold 
dark matter 

2.1. Physical Picture and Mathem,atical Formulation 

Assuming the linear growth of an overdense spherical top- 
hat perturbation, Gott (1975) and Bertschinger (1985) 
have derived that a bound object can grow, as a result 
of matter infall from the pressureless Hubble fiow, with 
the time dependence: M oc t^^^ and oc t^/^, where 
M is the total mass within the radius r„i of the growing 
bound object and t is the evolution time. Such a scaling 
relation depends only weakly on the decelerating param- 
eter go (for 0.02 < go < 2), and these relations become 
exact when go = 1/2. A slightly different scaling relation 
has also been predicted by Press and Schechter (1974) in 
their quest for the galaxy mass function. They found that 
M <x t^^^ and oc A feature common to both 

scaling relations is that the characteristic length scale Vm 
is approximately a linear function of time. This behavior 
will be used to fix the temporal scaling of velocity as i° in 
our self-similar model. 

On the other hand, the crossing of mass shells near 
the core can trigger particle orbit instabilities that yield 
violent relaxation (Lynden-Bell 1967). The particle orbits 
may be further randomized by the interactions with the 
crowded discrete galaxies in the cluster, making the dark 
matter particles further mixed. We therefore approximate 
the dark matter within the cluster to have an isotropic 
and also spatially uniform velocity dispersion a^,- Together 
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with the temporal scaling discussed in the last paragraph, 
one may consider Ud to be a constant in space and in 
time within the cluster. (For readers who have a different 
view of how the dark matter particles should behave in 
the cluster, they may simply regard the constant Od as a 
working hypothesis adopted to fix the self-similar scaling 
relations.) Outside the cluster, the cold dark matter de- 
tached from the expanding background Hubble flow falls 
radially inwards. The boundary of the hot cluster and cold 
background is defined by a heat conduction front, across 
which both mass flux and momentum flux are continuous. 

Thus, our picture of the self-similar matter infall can 
be summarized as follows. As the cold dark matter falls to- 
wards the massive compact core, a volume of finite radius 
surrounding the compact core then undergoes violent re- 
laxation and forms an isothermal halo, which gives rise to 
the cluster potential. A small fraction of the infall particles 
can however sink into the compact core and therefore the 
peak of the central density can grow steadily (Gott 1975). 
But the majority of infall particles are deflected back to 
the halo. The deflected particles find the halo to become 
more massive than they previously experienced due to the 
continuing infall, and hence the particles can travel only 
up to a maximum radius, which defines the location of the 
conduction front, i.e., the boundary of the hot cluster. As 
these processes continue, both compact core and halo can 
grow self-consistently. 

To put this picture into a quantitative perspective, we 
adopt a spherically symmetric infall model. Since the par- 
ticle orbits are assumed to be isotropic in the cluster bulk, 
the anisotropic stress force is unimportant and it suffices 
to consider only the radial pressure force in the momentum 
balance. The unique self-similar scaling with a constant (Td 
requires the mass density of dark matter particles p{r.t), 
the flow speed u{r,t) and the total mass M{r,t) enclosed 
within radius r to be scaled as: 
a{x) 



[{x-vf-l] 



P{r,t) 



u{r, t) = adv{x), M{r, t) 



G 



■m{x),{l) 
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where t = t and x (= r/adr) is the distance in a comov- 
ing frame. The continuity equation in the comoving frame 
then becomes 



da 



0. 



(2) 



— x^(2a + X—) + — (x^va) 
ax dx 

This equation can be rewritten as d{x'^{x—v)a)/dx = x^a. 
Comparing this equation with the relation between a and 
m, i.e., dm/dx = x^a, wc find that the mass m can be ex- 
pressed as m(.T) = x'^{x — v)a. This expression of m{x) is 
substituted into the gravitational force to eliminate m{x) 
in the momentum equation. Re-arrangement of the mo- 
mentum equation 



, dv 
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and eq.(2) yields 
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Note that eqs.(4) and (5) have been derived in the past 
for describing the inside-out collapse of a singular isother- 
mal molecular cloud (Shu 1977). However, our boimdary 
conditions are different from those for the molecular cloud, 
in that the hot cluster is fed from outside by the expand- 
ing cold dark matter, and inside, the coUisionless matter 
lands onto the compact core at a vanishing speed so that 
the particles can be captured. By contrast, the collapsing 
molecular cloud is fed from outside by a static isothermal 
gas sphere; inside, the collapsing gas lands on the pro- 
tostar violently with a Mach number much greater than 
unity, and it relics eventually on the radiative cooling for 
the infall matter to be captured by the prostar. 

Eqs.(4) and (5) can be simplified when the dark matter 
is cold, by setting the pressure term on the right of eq.(3) 
to zero. It yields 



1 da 
{x-v}- — 
a ax 



dv 
dx 

2. 

— [X — V 

X 



a 



(6) 
(7) 



These equations are valid for the background dark matter 
particles before falling within the cluster. Note that both 
eqs.(4) and (5), and eqs.(6) and (7) allow for the Hubble- 
flow solution of a flat universe where a = 2/3 and v = 
2x/3. 

Finally, the most relevant regime consistent with the 
current observations corresponds to a final saturated and 
static configuration when the matter infall ended before 
the present. That is, most of the available dark matter 
particles have collapsed into the cluster. Under these cir- 
cumstances, such a saturated system is completely deter- 
mined by the static forces in eq.(3): 



1 da 
a dx 



m 



(8) 



Here, wc have kept the same dimensionless parameter- 
ization although the dark matter density and mass are 
now time-independent. The additional factor r in this di- 
mensionless parameterization should be interpreted as the 
duration of cluster formation. 

2.2. Infall boundary conditions 

Our inner boundary condition of eqs.(4) and (5) is so cho- 
sen that a growing massive compact core is present. For 
the collapsed particles to settle into the compact core, we 
demand that the cluster bulk flow should have a vanishing 
infall speed at the core. Quantitative details as to how the 
matter falls into the compact core will be further elabo- 
rated below. 

Far from the cluster, we demand the pressure-free dark 
matter solutions to obey eqs.(6) and (7), and they must 
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either asymptotically match to a vanishing-density back- 
ground when ^l << 1, or asymptotically match to the 
Hubble flow (a = 2/3, u = 2x/3) when Q = 1. These 
pressure-free, outer infall solutions are then matched to 
the interior subsonic cluster solutions of eqs.(4) and (5) 
through a co-moving conduction front. Technically, we 
must adjust both the location and strength of the con- 
duction front to determine a unique solution satisfying 
both inner and outer boundary conditions. 

2.3. Solution in an open universe (fl <<1) or the 

ongoing infall case 

The background density of Hubble flow in this regime as- 
sumes a value a « 2/3, and our outer infall solution 
should asymptotically decrease to match this low-density 
background far from the cluster. Practically, this bound- 
ary condition demands the asymptotic outer infall solution 
to vanish at large distances. 

We shall stress that the solution to the open-universe 
case can also represent the situation where the self-similar 
infall began only recently in a flat universe, in which case 
T is redeflned as t — ti, with ti being the starting time of 
infall and r much less than the cosmic age t. As will be 
shown in the following analysis, the conduction front has a 
density p ~ I/AttGt^, much greater than the background 
density, e.g., the critical density Pc{= l/QwGt'^). There- 
fore the outer cold infall solution should asymptotically 
match the vanishingly small background density at large 
distances. Such a boundary condition is exactly the same 
as the open-universe case. We shall briefly return to this 
case in section 3.2 for a comparison with observations. 

The outer pressure-free solution is constructed from 
a remote location far from the cluster by integrating 
eqs.(6) and (7) inward. To ensure that the solution be- 
yond the starting distance of integration indeed asymptot- 
ically matches onto the background flow, we perform an 
asymptotic analysis for the large-distance solution. The 
large-distance solution of eqs.(6) and (7) is a = cq/x'^ and 
V = —cq/x. That is, the density a and velocity v must 
be related in this specific way to warrant that the density 
asymptotically decreases to a sufficiently low value and 
matches onto the low background density. However, since 
eqs.(6) and (7) are actually nonlinear equations, the am- 
plitude Co must also be constrained in such a way that the 
outer infall solution be matched to the interior hot cluster 
solution through a conduction front. In this regard, the 
amplitude cq can be taken as a nonlinear eigenvalue im- 
posed at the outer boundary. Moreover, the interior clus- 
ter solution is demanded to possess a diminishing velocity 
when the flow enters the compact core, so that the infall 
matter can be captured by the growing compact core. This 
requires a unique location for the conduction front as the 
second nonlinear eigenvalue. 

To meet these non-trivial conditions, we find it useful 
technically to search for an approximate solution of eqs.(4) 



and (5), which makes a smooth transition from a subsonic 
infall in the bulk to a supersonic infall into the core. This 
is because the trans-sonic solution is at the bifurcation 
boimdary in the solution space (Shu 1977), and our de- 
sired solution is in the neighborhood of this boundary. 
The desired solution can then be constructed relatively 
easily after the approximate solution is found. Conditions 
for a flow to smoothly pass through the sonic speed are 
as =2/xs and Vs = Xs — l, where the subscript s denotes 
the sonic point which is a critical point of eqs.(4) and (5). 

By adjusting both the amplitude cq of the asymptotic 
background solution and the location Xs of the sonic crit- 
ical point, one may uniquely determine the location of 
and density jump at the conduction front, across which 
the mass flux and momentum flux in the co-moving frame 
are continuous. Our numerical search shows that the ap- 
proximate solution has cq « 1.22 and Xg « 10""^, and the 
location of the conduction front at Xcond ~ 2/3. Upstream 
of the conduction front, the dark matter density a « 1.6 
and infall speed v « —1.45, whereas downstream of the 
front, a w 3.5 and v w —0.32. This conduction front has 
a moderate compression ratio « 2.2. 

An eigenvalue search in the neighborhood of this trans- 
sonic solution permits us to obtain the desired solution 
easily, which turns out to differ from the trans-sonic solu- 
tion only near the compact core. As previously stated, the 
transonic solution is the bifurcating solution, and it has an 
infall speed increasing inward moderately. On one side of 
this bifurcation boundary, the infall speed increases dras- 
tically inward, but on the other side the infall speed turns 
around to become vanishingly small toward the core. (See 
Fig. 1.) (These situations are not unlike the behaviors of 
the scaling factor a{t) in a flat universe, an open universe 
and a closed universe.) Clearly, the latter is what we look 
for. As the desired solution and trans-sonic solution are al- 
most identical in most part of the cluster, their conditions 
at the conduction front are practically the same. 

A caveat for our results is that the coUisionless conduc- 
tion front is likely to have a transition layer of finite thick- 
ness. Hence the results given above can only be regarded 
as a sophisticated estimate for modeling the boundary of 
the cluster, and the discontinuity obtained here should not 
be taken too literally. 

The derived density profile of dark matter fiow can be 
fit nicely by an analytical form: 

{(0 13^-|-x^)^/^ ^ ^ ' "^^ X <C Xcond'i 

(9) 
Q 52^^+3;^ ' > X > Xcond- 

Two separate components for the interior cluster solu- 
tion are identified. The first component (oc e^-^*^^^/^) gives 
rise to a peaked density profile, in response to the gravity 
of massive compact core; the second component represents 
the usual Bonner-Ebert state of isothermal sphere, which 
has a secondary core of size Xc ~ 0.13, and the density out- 
side the outer core decreases outward first as x~^ and sub- 
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sequently as The conduction front is located roughly 
at the transition between the and profiles. The 
matter located in the region from x ~ 0.02 to a; 0.4 can 
be well approximated by a "static sphere" , where the flow 
speed is roughly zero. Near the conduction front the infall 
speed sharply rises, and the finite flow speed immediately 
downstream of the conduction front serves only to flatten 
the density profile slightly. At about x = 0.02, the mass 
within this sphere is twice as much as that of the compact 
core; that is, the gravity arising from the massive com- 
pact core begins to be overwhelmed by the self-gravity of 
the matter outside the compact core. Figs. 1(a), (b) and 
(c) depict the normalized density, velocity and mass pro- 
files, respectively. Note that most part of the hot cluster 
is almost static, except near the conduction front and im- 
mediately outside the compact core. 



2.4- Infall into the compact core 



The compact core is a virialized bound object, where the 



lowing analytical expression is found to fit the outer cold 
infall solution: 



density rises inwards more stc!c;ply than 



As shown 



in Fig.(l) when the dark matter particles fall within x = 
0.02, they experience solely the core gravity and are accel- 
erated inward to reach a maximum speed at x ~ 8 X 10"'*. 
Inside x ~ 8 x 10"'', the infall suddenly gets retarded by 
the core pressure and its speed decreases from this point 
on and becomes vanishingly small at the core center. 

To pin down where the true core boundary is, we ex- 
amine the binding energy. Inside the core, the fiow must 
satisfy the binding condition, \P-E.\ > 2\K.E.\, so that no 
particle may escape outside the radius x once enters. The 
specific potential energy P.E. at the radius x is —m/x, and 
the specific kinetic energy K.E. includes the specific ther- 
mal energy 3/2 and the comoving flow energy {v — x)^/2. 
Our solution shows that only inside x = 5 x 10~^ can the 
binding condition be satisfled. The size of the compact 
core is thus set to be a;6 = 5 x 10"^. 



2.5. Solution of early infall in a flat universe 

The asymptotic solution far away from the cluster behaves 
as V = (2a;/3) — ci/x and a = (2/3) -|- ci/x^, when the 
deviation from the Hubble flow is small. We may employ 
the technique used in the previous case to obtain the solu- 
tion in this regime. Surprisingly, the trans-sonic solution 
yields a conduction front at almost the same location and 
with a comparable jump as those of the previous case. The 
nonlinear eigenvalue Xg for the trans-sonic solution is also 
the same. In addition, the nonlinear eigenvalue Ci also as- 
sumes a value equal to ci = 1.22. As a result, we flnd that 
the desired solution interior of the cluster is nearly the 
same as that of the previous case. The outer cold infall 
solution is, however, different from that in eq.(9); the fol- 



a{x) 



0.91 2 



10> X > X. 



cond- 



(10) 



2.6. Saturated phase 



Taking the typical dark-matter velocity dispersion, a^. « 
1000 km s~^, and the typical size of an Abell cluster, 
3 Mpc, it follows that the typical time scale for cluster 
growth is on the order of 10^ years, a short time compared 
with the age of the universe. Moreover, no strong cosmic 
evolution for clusters since z ^ 0.8 has been detected by 
recent X-ray observations (Rosati, et al. 1988; Vikhlinin 
et al. 1998). It is, therefore, reasonable to assume that 
most clusters have reached saturated states by z « 0.8. 
The comparison of observations with the present model is 
to be made primarily for the saturated state, where the 
clusters arc virialized. 

Because the central compact core is strongly gravita- 
tionally bound, the core should remain intact in the sat- 
urated phase, and the core size and mass are frozen to 
their values immediately before the infall stopped. Note 
that the extra factor r that makes eq.(7) dimensionless 
can be interpreted as the cosmic age at which the infall 
ended for an early formed cluster. Therefore, we may con- 
veniently choose the same values of compact core mass 
(t77 = 1.6 X lO^'') and size {xi, = 5 x 10^"') as those of the 
infall phase for the static cluster solution. 

The integration of eqs.(2) and (7) can be carried out 
from the compact core outward to a sufficiently large dis- 
tance to obtain the static cluster solution. The solution 
construction is straightforward once the core size and mass 
are specified. The resulting density profile of dark matter 
is plotted in Fig. 1(d) and can also be fit by an analytical 
form: 



aix) = 1.68 e°-0°i5/^ 



(0.122 +a;2)3/2" 



(11) 



The density near the core also behaves as cx e" °"i''^/^", and 
again has a secondary outer core Xc = 0.12, beyond which 
the dark matter decreases outward also from a oc x~^ pro- 
file to oc .x~2 profile at large distances. The corresponding 
profiles are depicted in Fig.l(c & d). 

In sum, for the solutions of both ongoing- infall (open 
universe) and early- infall (flat universe) cases, the inte- 
rior hot clusters are found to be almost identical. This re- 
sult to some degree reinforces the assertion of Gott (1975) 
that the self-similar evolution of bound objects depends 
weakly on the deceleration parameter. In the saturated 
phase, we have assumed that the compact core is frozen 
to the conflguration immediately before the infall ended, 
and obtained a saturated cluster configuration that de- 
parts only slightly from those of the infall phase. We also 
find that the plausible value for the compact core radius 
approximately equals Xf, « 5 x 10~^ and the core mass. 
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Fig. 1. (a)Diniensionless density profile, (b)flow speed and (c) cluster mass obtained from the recent-infall, self-similar 

solutions under the initial and boimdary conditions of a mass concentration mo = 1.6x10"^ within x = r/adT = 0.0005. 
Also shown are the density (d) and mass (c) variations for a saturated cluster (dashed lines) expected under the same 
inner boundary condition. The dotted lines in (d) are the analytic fits of the density profiles by eqs.(9) and (10), 
respectively. 



1.6 X 10~^. Taking a typical cluster for which x = 1 corre- 
sponds roughly to 1 Mpc and a typical velocity dispersion 
~1000 km s~^, we estimate Xb ~ 0.5 kpc and the compact 
core mass about 3 x 10^ ^Mq. 

3. Comparison with observations 

3.1. Rationale and strategy 

In this section, the predicted dark-matter distribution in 

the galaxy cluster will be tested against the observed re- 
sults of gravitational lensing, galaxy distribution, gas dis- 
tribution and baryon fractions. The rationale behind such 
a comparison is that if the theoretical prediction can well 
approximate the true dark-matter distribution, a single 
gravitational potential given by the predicted dark-matter 
distribution ought to be able to produce consistent pro- 
files for these different observations. Our strategy is there- 
fore first to use the gravitational lensing data to calibrate 
the length scale and depth of the gravitational potential, 
which are related to the only two parameters, r and aj,, 



in the theoretical model. Once the potential is fixed by 
the lensing data, the baryons will be treated as test parti- 
cles, which should be so distributed as to be in dynamical 
equilibrium with the potential. 

However, the construction of galaxy and gas distribu- 
tion also requires an additional knowledge of their "equa- 
tions of state" . For galaxies, we adopt the average profile 
of the observed galaxy velocity dispersion and for g ases, a 
polytropic equation of state. Though it can be relatively 
easy to find the best galaxy and gas parameters for a the- 
oretical model to compare well with the observations, we 
would like to go one step beyond, and to constrain the 
choice of these empirical parameters based on a physical 
principle. That is, we consider the dark matter particles, 
galaxies and gases to have already been grossly virialized, 
and they should share a common averaged velocity dis- 
persion since they are all within a common gravitational 
potential. It is under this additional constraint that the 
comparisons will be made. 
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3.2. Dark matter distribution: comparison with the 
gravitational lensing results 

Gravitational lensing appears to be a powerful and unique 
tool at present for probing the total matter distribution of 
galaxy clusters. It provides the gravitating masses of clus- 
ters regardless of the matter compositions and dynamical 
states. So far, the strongly and weakly distorted images 
of background galaxies have been detected in about 50 
clusters at intermediate redshifts ranging from .2 « 0.1 to 
z « 0.8. These clusters constitute a good sample of mat- 
ter distributions on scales r ~ 10 kpc - ~ 1 Mpc, both 
individually and statistically for testing various mass den- 
sity profiles advocated by theoretical and observational 
studios. A direct comparison between om' solution eq.(ll) 
and the gravitational lensing results is straightforward, 
provided that the evolutionary time r and the velocity 
dispersion ad of dark matter particles arc known. On the 
other hand, such a comparison will allow us to fix these 
two free parameters, if the self-similar solution predicts 
a consistent mass profile with what is derived from the 
gravitational lensing. Due to the cosmic variance, these 
two parameters can only be fixed sensibly by the aver- 
aged mass profile of the lensing results. 

Before we fix these parameters by the averaged lensing 
results, a comparison with two individual clusters is pre- 
sented to illustrate how well the theoretical profile may 
agree with the observed profiles individually, and not just 
statistically. Fig. 2 shows the projected radial mass distri- 
butions of two distant clusters at z « 0.8, MSI 137-1-66 and 
RXJ1716+67, obtained recently by employing the mass re- 
construction technique for the detected weak distortions 
of distant background galax;ies (Clowe et al. 1998). Also 
plotted are the 2-D mass profiles from eq.(ll) for ((7^, 
(t^t)=(850 km s-\ 1.5 Mpc) and {ad, crdT) = (1000 km 
s~^, 0.6 Mpc), respectively. (From now on, the quantity r 
should, unless otherwise indicated, be interpreted as the 
cluster formation time scale.) It appears that the self- 
similar solution does yield an acceptable fit to the results 
of MS1137-I-66 and RXJ1716+67. Moreover, the required 
values of velocity dispersion, ad ~ 1000 km s~^, for dark 
matter are also in good agreement with the observed X- 
ray luminosity of the two clusters, Lx ~ (6 — 8) x 10**^ 
erg s~^ in the energy band 0.3-3.5 keV. Meanwhile, the 
scale length adT ~ 1 Mpc indicates a comfortable cluster 
formation time of r ^ lO'' yrs. 

This apparent success encourages us to conduct a com- 
parison between the overall radial mass distribution of 
clusters revealed by the strong/weak lensing and that pre- 
dicted by our analytical solutions. Note that in most cases, 
only the rich and X-ray luminous clusters of galaxies at in- 
termediate redshifts are capable of gravitationally distort- 
ing the background galaxies. Namely, the (strong) lens- 
ing clusters alone are a fair sample of massive clusters at 
z « 0.3 in the universe. We display in Fig. 3 the mean pro- 
jected cluster mass miens within r, derived from a total of 




0.5 1 1.5 



r (Mpc) 

Fig. 2. A comparison between the projected cluster 
masses derived from the weakly distorted images of 
background galaxies behind clusters MS1137-I-66 and 
RXJ1716-I-67 and those from our saturation model. Data 
sets of the lensing cluster masses are from Clowe et 
al. (1988). Our best fitting yields {ad,adT) = (850 km 
s-1,1.5 Mpc) and {ad,adT) = (1000 km s"i,0.6 Mpc) for 
MS1137H-66 and RXJ1716-I-67, respectively. 

49 arcs/arclets in 38 clusters and 134 weak lensing mea- 
surements made among 24 clusters (Wu et al. 1998 and ref- 
erence therein), in which we have not included the result of 
A2163. While there may exist a discrepancy between the 
strong and weak lensing measured cluster masses, the clus- 
ter masses obtained from three different lensing phenom- 
ena, including the deficit of red galaxy population detected 
behind A1689 due to the lensing magnification (Taylor et 
al. 1998), are essentially consistent within 2a error bars. 
The radial mass variation of the lensing clusters from our 
self-similar solution with ad = 1300 km s~^ and adT = 1 
Mpc is also demonstrated in Fig. 3. The good fitness of our 
prediction to the lensing result is clearly seen. Since the 
velocity dispersions of optical galaxies in lensing clusters 
are in the range of 600 - 2000 km ~^ with an average value 
of 1200 km s~^ (Wu et al. 1998), the required mean veloc- 
ity dispersion of dark matter particles ad = 1300 km s~^ 
is more than acceptable. Consequently, the mean forma- 
tion time of the lensing and massive clusters turns out to 
be r = 7.5 x 10^ yrs. This timescale seems to be relatively 
short though it is not impossible for the very rich clusters 
(Gunn & Gott 1972). 

With the same velocity dispersion and scale length 
parameters, our theoretical results produce equally good 
fits for the mass profiles given by both the static solu- 
tion of eq.(ll) and the collapsing solution of eq.(9). The 
latter should now be interpreted as the solution of a grow- 
ing cluster which began its formation only recently. Also 
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Fig. 3. An overall radial variation of the projected clus- 
ter masses revealed by the gravitational lensing methods: 
the strongly (arcs/arclets) and weakly distorted images of 
background galaxies (for references sec Wu et al. 1998) as 
well as the deficit of red galaxies behind A1689 (Taylor et 
al. 1998). The results predicted by the recent-infall solu- 
tion (dotted line) and the saturated configuration (solid 
line) are shown for (Jd = 1300 km s~^ and g^t = 1 Mpc. 

shown in Fig. 3 is that the difference between the mass 
profiles is only minor, leading to an inability to observa- 
tionally distinguish the two models simply from their mass 
profiles. To differentiate these two models, it is helpful to 
recall some established observational evidences of galaxy 
clusters. It was shown more than a decade ago that op- 
tical counts of clusters did not exhibit any evolutionary 
tendencies for redshift out to at least z « 0.5 (Gunn, 
Hocssel, & Okc 1986). The same conclusion holds true 
for the X-ray selected clusters since redshift of as high as 
z 0.8 (Rosati et al. 1998 and references therein). Fur- 
thermore, no significant differences in the dynamical prop- 
erties between low-redshift and high-redshift clusters, such 
as the X-ray luminosity-temperature relationship, the ve- 
locity dispersion-temperature relationship (Mushotzky & 
Scharf 1997; Wu et al. 1999), and the baryon fraction (see 
Fig. 7), have been detected. In particular, the distribution 
of core radii of the intracluster gas in nearby clusters is 
identical to that of distant clusters {z > 0.4) (Vikhlinin et 
al. 1998). This last point can be regarded as the most con- 
vincing evidence for the already "settled" configuration 
of cluster matters, but sharply discords with an evolv- 
ing self-similar model, which predicts the core radius to 
be much smaller in the past (Kaiser 1986). Suggested by 
these evidences can be that significant matter infall pro- 
cesses should have already ended before z 0.4 so that 
the global properties of galaxy clusters have remained ap- 
proximately unchanged since an even earlier epoch. 



3.3. Galaxy distribution 

The rationale of the following exercise goes as follows. 
From the massive-cluster samples measured by gravita- 
tional lensing, we have statistically determined the mean 
velocity dispersion ad{~ 1300 km s"'^) and the mean 
length scale aar (~ 1 Mpc), thereby fixed the average 
mass and length scales of our solutions for massive clus- 
ters. This predicted average sohition hence has no ad- 
justable parameter. So, we may employ the Jeans equation 
to construct an expected average galaxy distribution and 
compare it with the observed average galaxy distribution 
within the massive clusters. 

So far we have not yet included the baryonic matter 
(galaxies and gas) in our self-similar solutions of dark mat- 
ter particles in clusters. Although an exact self-similar so- 
lution with a proper combination of the collisionless dark 
matter and the coUisional gas infall can be constructed, we 
would rather take a less vigorous approach to this problem 
by treating the baryons as test particles. For galaxy pop- 
ulations characterized by their radial velocity dispersion 
o'gaii'i^) and number density Ugaiir), the Jeans equation 
under the light-traces-mass hypothesis reads 

1 dialiUgai) ^ 2al^^^_ GM ^ 
Ugai dr r ' 

where M is the total cluster mass (we use dark matter in- 
stead) enclosed within radius r, and A denotes the veloc- 
ity anisotropy parameter. Since the above equation corre- 
sponds to a purely hydrostatic equilibrium, we employ the 
static solution, cq.(ll), to evaluate M{r). (The collapsing 
solution, eq.(9), has also been tested, and essentially the 
same M(r) is obtained.) Once the velocity profile <Tgai{r) 
and A are specified, we can easily work out the galaxy 
number density profile according to eq.(12). 

We compare our prediction with the average surface 
number density of galaxies over 14 distant clusters mea- 
sured by the Canadian Network for Observational Cosmol- 
ogy (CNOC) (Carlberg, Yee, & ElHngson 1997). Actually, 
the majority of the 14 clusters arc those massive lensing 
clusters discussed in the above subsection. We adopt the 
simple form of the radial velocity dispersion suggested by 
CNOC, <yg^i {r) = (Tg„((0)/(l-|-r/r2oo) and the assumption 
for a zero anisotropy parameter, A = 0. Here, agai{0) is 
the central velocity dispersion and r2oo denotes the radius 
within which the mean cluster mass density is 200 times 
the critical mass density of the universe {Qq = !)• The 
mean value of r2oo over the 14 CNOC clusters has been 
found to be r2oo = 2.5 Mpc. 

With the predicted gravitational potential of the un- 
derlying dark matter and the given profile of observed 
velocity dispersion of galaxies, one can solve for the aver- 
age galaxy distribution using the Jeans equation, eq.(12). 
In the following comparison, we will instead take an al- 
ternative route by first seeking the best fit parameters 
and then examine how well they may agree with those 
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Fig. 4. Surface number density profile of galaxies over an 
ensemble of 14 distant clusters from CNOC (Carlberg, 
Yee, & Ellingson 1997). Solid line is our expected re- 
sult of the saturation model for (Jd/(Jgai{0) = 0.78 and 
i'2Qol<ydT = 0.4 under the assumption that galaxies trace 
the gravitational potential of dark matter particles in clus- 
ters. The conventional King model fit is illustrated by the 
dotted line. 

fixed by the lensing results. Fig. 4 shows the CNOC galaxy 

surface number density Sg^; and the predicted galaxy 
surface number density from our model for the best fit 
fd/fga/(0) = 0.78 and cFd.T/r2oo = 0.4. When combined 
with r2oo = 2.5 Mpc, this choice of length scale (TdT{= 
0.4r2oo = 1 Mpc) is what has previously been determined 
by the gravitational lensing. The good fit of the predicted 
length scale with the observed length scale is remarkable. 

The magnitude of the predicted surface number den- 
sity depends sensitively on the central velocity bias 
(Td/(Tga;(0). Its best fit value being less than unity may sim- 
ply reflect the fact that the central galaxy velocity agai (0) 
is higher than the average value agai- To pin down the 
ratio of the averaged specific kinetic energy of galaxies to 
that of dark matter particles for this value of <Jgai{0), we 
calculate 



= cr, 




T,gal{x)xdx 



(13) 



A straightforward calculation shows that Odjogai ~ 1- 
That is, the best-fit central velocity bias (Td/o-gaKO) ren- 
ders the average specific kinetic energies of both dark mat- 
ter particles and galaxies to be approximately equal, as re- 
quired by our principle set forth to constrain the empirical 
parameters. The good agreement is more than surprising 
for such a simple theoretical model. Actually, it has al- 
ready been seen in the above subsection that the velocity 
dispersion of dark matter particles (cr^ = 1300 km s~^). 



required to account for the projected mass distribution of 
clusters derived from gravitational lensing, is roughly the 
same as the reported mean value for the CNOC galaxies 
(o'gaj = 1200 km s"-'^)). This is also in agreement with the 
previous results given by a statistical comparison of cluster 
mass determinations from gravitational lensing and veloc- 
ity dispersion of cluster galaxies as the tracer of cluster 
potential (Wu & Fang 1997; Wu et al. 1998). 

3.^. Gas distribution and (3 model 

Similarly, the equation of hydrostatic equilibrium can also 

be applicable to the hot X-ray emitting gas of number den- 
sity ngas{r) in clusters (Cavaliere & Fusco-Femiano 1976): 



1 dial 



s) 



GM 



n, 



gas 



dr 



(14) 



in which (7^^^ = kT/iinip is the "equivalent velocity 
dispersion" of gas with fxrUp being the mean particle 
mass. Assuming an equation of state for the gas T = 
To [n gas /n gas (O)]''^ ^ , we can write out the gas number den- 
sity profile in terms of eq.(14): 



Tin 



1- 



7-1 

7 



Jr. ^. 



GM 



1 

gas 



(0)r 



-dr 



1/(7-1) 



(15) 



where rf, is the physical length corresponding to xj-,. Moti- 
vated by the fact that 'Ogai ~ cr^, the arbitrariness of the 
polytropic index 7 and (7gas{Q) can be fixed by a physical 
requirement that 'a gas also equals Ud 
average equivalent velocity dispersion of gas: 



where Ugas is the 



{x)S{x)xdx 



S{x)xdx 



(16) 



The observed X-ray surface brightness is related to the 
temperature and number density of electrons in clusters 
as 



T'/\x)nl 



xdx 



~2' 



(17) 



provided that the radiation mechanism is the thermal 
Bremsstrahlung with a Gaunt factor weakly dependent 
on T. 

Before proceeding to a discussion on the general prop- 
erties of the theoretically predicted X-ray surface bright- 
ness profile, we first present a comparison of the predicted 
profile with the observed X-ray surface brightness of an 
individual cluster for an illustration of the agreement that 
also holds individually. Figure 5 shows a typical exam- 
ple of the i?05^r observed S^{r) for cluster C10016-hl6 
[z = 0.55), together with our expected Sx{x) from eqs.(15) 
and (17). The latter are obtained by properly choosing the 
length scale parameter cjrfT and the central biasing pa- 
rameter <Td/agas{0) or equivalently the polytropic index 
7. This procedure is somewhat equivalent to the conven- 
tional P model fit by properly determining the X-ray core 
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Fig. 5. An examples of fit of tlic theoretically predicted 
X-ray surface brightness (solid lines) to the observation 
of C10016+16 (Neumann & Bohringer 1997). The best /3 
model fit is also displayed (dotted lines) for comparison. A 
background surface brightness of 5.5 x 10~^ s~^ arcsec"^ 
is assumed. 



radius and the power index. Yet, unlike the /3 model which 
is purely an empirical formula, our predicted X-ray surface 
brightness profile is given by the hydrostatic equilibrium 
of the intracluster gas with the underlying gravitational 
potential of the cluster dark matter particles. The physics 
behind our prediction about the gas distribution is well 
understood. 

Fig. 6 shows the X-ray surface brightness profiles of 
clusters obtained from eqs.(15) and (17) for a set of poly- 
tropic indices 7 =1.05, 1.1, 1.2 and 1.3, in which wc 
have imposed a constraint on Od/ <ygas{^) by requiring 
^d/'^gas = 1) which Uniquely defines a central biasing pa- 
rameter (Jd/ <^gas{'^) for each value of 7. Wc now superpose 
the conventional /3 models on the theoretically expected 
Sx- Surprisingly, all the derived profiles can be well 
approximated by the /3 models except near the compact 
cores, where brightness cusps may appear in some cases. 
The allowed /3 parameters are dependent on the adopted 
7, with a smaller (3 for a large 7. Additionally, our fit- 
ted (3 parameters are consistent with the observed (3 pa- 
rameters in the range of /3/it ~ 0.6 - 1.2. If the average 
value of (3 fit is indeed aroimd 2/3 as indicated by most 
current observations, we may conclude that a large pro- 
portion of galaxy clusters should have a mean polytropic 
index of 7 « 1.2-1.3. This point has recently been verified 
by the X-ray measurements of temperature structures of 
30 nearby clusters (Markevitch et al. 1998), all of which 
show similar temperature profiles with a polytropic index 
of 7 = 1.2-1.3. 



On the other hand, the above fits for various 7's 
yield that the core radius of the gas profile varies from 
Tc/f^dT ~ 0.16 to 0.30. Of these, the choice for 7 = 1.2 
that yields TcjodT ~ 0.20 is again of particular interest. 
Since if the mean length scale OdT is approximately ~ 1 
Mpc for rich clusters, as was shown in the above two sub- 
sections, this choice implies a mean core radius of ~ 0.2 
Mpc in Sx'-, this core radius is comparable to the median 
values (0.21-0.24 Mpc) of the X-ray core radius distri- 
butions of the ROSAT detected clusters (Vikhlinin et al. 
1998). Taking these results as a whole, we have provided 
a consistent scenario for the physical origin of the conven- 
tional f3 model. 

In sum, our fit suggests that the gas components 
of clusters have an overall polytropic index close to 
7 w 1.2 and a mean (3 parameter 2/3. It 

is worth stressing again that by requiring the ratio of 
mean specific kinetic energy in galaxy to that in gas 
to be unity, f3spec = wl^i/{kTgas/ ump) = ^l^i/Trlas = 
{f^dlogasf lipdjogaif' « 1, the two /? valucs, (3fn and 
(3 spec, turn out to agree nicely with the observationally 
determined ones (Wu, Fang & Xu 1998). The two /3's are 
not the same, and the so-called "/? discrepancy" (Bahcall 
& Lubin 1994; references therein) appears to be a natu- 
ral consequence of our model. We have also noticed that 
a similar effort was recently made by Makino, Sasaki & 
Suto (1998) based on the universal dark matter profile. 
However, while their derived X-ray surface brightness is 
well represented by the /3 model in shape, their predicted 
core radius is a factor of (3-10) smaller than the actually 
observed values. 

3. 5. Baryon fraction 

Having known the dark matter and gas distributions, we 
can easily provide the baryon (gas) fractions of clusters 
fb = Mgas{r)/Mtoiai{r). where Mgas and Mm al arc the 
gas and total masses enclosed within the radius r, respec- 
tively. Of course, the absolute value of fb is ultimately 
related to the amplitude of gas density profile, ngaj;(0), 
which has remained unknown thus far. Although we shall 
be able to figure out Ugas (0) using the observationally de- 
termined baryon fractions, a particular attention is paid 
to the issue of how the baryon fractions vary with scales. 
Fig. 7 demonstrates the cluster baryon fractions versus 
radii obtained through the deprojcction analysis of Ein- 
stein Observatory data for 207 clusters (White, Jones & 
Forman 1997). Basically, the baryon fraction is an increas- 
ing function with radius and tends towards an asymptotic 
value of ~ 20%-30% at large radius. Such a variation 
has been noticed previously by a number of authors (e.g. 
White & Fabian 1995; David 1997; White et al. 1997) but 
not yet been satisfactorily accounted for in the framework 
of any prevailing models for cluster formation. 

Adopting a mean polytropic index of 7 = 1.2 for the 
X-ray gas, we attempt to reproduce the observed radial 
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Fig. 6. X-ray surface brightness profiles predicted for different polytropic indices 7 = 1.05 1.1, 1.2 and 1.3 in tlie 

equation of state T oc n^~^ (solid lines). The requirement that the mean kinetic energy of dark matter particles is 
equal to that of gas, cr^ = {kT / jimp) has been employed. Also plotted are the conventional /? gas profiles (dotted 
lines). 



variation of baryon fractions with our model. Since the 
length scale of our model has already been fixed, the fit- 
ting can only be conducted in terms of the absolute value 
of baryon fractions. The best-fit curve to is illustrated in 
Fig. 7. The significant rise of /& with radius from r ~ 0.1 
Mpc outward is due to the fact that in this radial interval 
the dark matter profile goes as (Fig.l) whereas the 
gas distribution varies as (Fig. 6c). However, it should 
be pointed out that the validity of our prediction about 
outside r ^ 5 Mpc begins to be questionable, though a 
decreasing tendency is seen to occur beyond r ^ 10 Mpc. 

Finally, we have noted that the majority of the data for 
the baryon fractions fb have so far been derived from the 
isothermal gas hypothesis, for which an emission-weighted 
temperature over a finite X-ray emitting surface is applied. 
Therefore, the comparison between our prediction with 
the measured baryon fractions holds true only statistically. 



4. Discussions and conclusions 

In the present work, we explore the possibility of building 
up two distinct mass cores in a single galaxy cluster as a 
result of matter infall from the background Hubble flow. 
Comparisons of the predicted saturated configuration with 
observations reveal that our model agrees quite well with 
the lensing, galaxy distribution and gas distribution. Our 
result is insensitive to the density parameter since the 
cluster builds up within a period much shorter than the 
age of the universe. The key assumptions of this model 
are (a) a pre-existing mass concentration, onto which the 
background matter falls, (b) the infall within the cluster 
is subsonic and the dark matter particles are largely viri- 
alized with ergodic orbits, in which situation the particle 
velocity dispersion is assumed to be a constant and (c) 
in the saturated phase the mass ratio of the two cores 
remains the same as that in the infall phase. 

Despite of its oversimplification, the self-similar grav- 
itational collapse of collisionless dark matter particles re- 
mains to be the only analytic approach to access the dy- 
namical process of cluster formation, which may comple- 
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Fig. 7. The binned radial variation of baryon (gas) frac- 
tions (/(,) of clusters. The data sets are from the Einstein 
Observatory (White, Jones, & Forman 1997; WJF) and 
from the measurements of the S-Z effects in four clusters 
(Myers et al. 1997). The high-redshift clusters {z > 0.1) 
in WJF are explicitly plotted for comparison. The allowed 
cosmological baryon fraction in the; standard Big Bang 
Nucleosynthesis (BEN) for a matter dominant flat uni- 
verse is also illustrated. The solid line is the theoretically 
predicted baryon fraction by our saturation model for gas 
with 7 = 1.2 and ad/a-gas{0) = 0.62 

ment our understanding of its physical causes in addition 
to the employment of large-scale numerical simulations. 
In this paper we have re-addresscd the self-similar infall 
model in a new way: Our initial and boundary conditions 
are so chosen that there is a pre-existent mass concen- 
tration in the central position of cluster to induce mass 
accretion, and the process of major dark matter infall is 
confined within a finite interval of ~ 10^ years. Conse- 
quently, a hydrostatic equilibrium is established by now 
between the gravitational potential of dark matter and 
the galaxies/gas in the late evolutionary stage of clus- 
ters (probably since rcdshift z ^ 0.8). In the scenario 
that the global properties of clusters remain not evolv- 
ing, we have successfully reproduced the distributions of 
total mass, galaxies and hot gas in clusters of galaxies 
revealed by gravitational lensing and optical/X-ray obser- 
vations. Apart from the very central massive core of ~ few 
lO^^M;;.) within the radius of 1 kpc (adopting ad = lO'' 
km s^^ and a^r = 1 Mpc), the dark matter, galaxies 
and gas essentially share a common "secondary" core at 
r ~ 0.15 — 0.2 Mpc, though their profiles can differ signif- 
icantly, depending on the radial variations of the velocity 
dispersion/ thermal speed of galaxies/gas. 

We also find that the bias of the cluster averaged ve- 
locity dispersion of galaxies against the velocity disper- 



sion of dark matter is indeed approximately imity, i.e., 
o'gai/o'd ~ 1, as expected from the virialization considera- 
tion that they are both within the same gravitational po- 
tential. In view of this result, it is further demanded that 
the specific dark matter kinetic energy and specific gas 
thermal energy to also be equal, i.e., {kT gas / tJ-mp) / cr'^^ w 1, 
thus yielding the /3spec parameter, which denotes the ratio 
of specific kinetic energy in galaxy to that in gas, to be also 
approximately unity (Bahcall & Lubin 1994; Girardi et al. 
1996). The resulting X-ray surface brightness is foimd to 
be well fit by the conventional /3-model, with, however, the 
parameter ^spec different from the (ifn parameter used in 
the /3-modcl fitting of the X-ray surface brightness. The 
mean value of (jfu turns out to be /3jjj « 2/3, in contrast 
to /Jspec ~ 1- Therefore, the (3 discrepancy, li spec / 1^ fit 
arises, in our model, as a natural consequence of the poly- 
tropic index of gas 7 7^ 1. With the mean value /3jjt 2/3, 
it yields that the mean value of 7 should be 7 w 1.2. 

Apparently, the present model has its own problems 
that remain to be resolved in the future. First of all, it 
is unclear at what epoch the central massive core should 
begin to form; this issue is relevant to the question as to 
when our model starts to be applicable. Second, although 
optical/X-ray observations prefer a non-evolutionary sce- 
nario of galaxy clusters within z ^ 0.8, we are unaware of 
how the background matters should terminate the infall 
onto clusters. It is possible that after the cluster mass 
reached some threshold value the intercluster matters 
could then be heated up by the intracluster tidal forces 
exerted by clusters, thereby preventing further matter in- 
fall. How this heating mechanism may proceed remains 
to be examined. Finally, the massive compact core is pre- 
dicted to be about 1 kpc in size. However, the current 
observations arc unable to set any robust constraints on 
the matter distribution inside the core region as small as 
r < 10 kpc. It requires future observations to confirm our 
prediction of compact cores at the cluster centers. 
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